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The atomic ionization processes induced by scattering of neutrinos play key roles in the experimental 
searches for a neutrino magnetic moment. Current experiments with reactor (anti)neutrinos employ germa- 
nium detectors having energy threshold comparable to typical binding energies of atomic electrons, which 
fact must be taken into account in the interpretation of the data. Our theoretical analysis shows that the 
so-called stepping approximation to the neutrino-impact ionization is well applicable for the lowest bound 
Coulomb states, and it becomes exact in the semiclassical limit. Numerical evidence is presented using the 
Thomas-Fermi model for the germanium atom. 



The neutrino magnetic moments (NMM) expected 
'in the Standard Model are very small and proportional 
to the neutrino masses [T]: /i^ « 3 x 10^^^/is(m,y/l eV) 
.with /is = e/2m being the electron Bohr magneton, 
'and m is the electron mass. Thus any larger value of 
/Zj/ can arise only from physics beyond the Standard 
Model (a recent review of this subject can be found in 
'Ref. [2]). Current direct experimental searches [3l|4l[5] 
|for a magnetic moment of the electron (anti)neutrinos 
'from reactors have lowered the upper limit on ij,,^ down 
.to < 3.2 X 10~^^/iB [5]. These ultra low background 
'experiments use germanium crystal detectors exposed 
to the neutrino flux from a reactor and search for scat- 
tering events by measuring the energy T deposited by 
the neutrino scattering in the detector. The sensitivity 
of such a search to NMM crucially depends on lower- 
ing the threshold for the energy transfer T, due to the 
enhancement of the magnetic scattering relative to the 
standard electroweak one at low T. Namely, the dif- 
ferential cross section da/dT is given by the incoherent 
sum of the magnetic and the standard cross section, and 
for the scattering on free electrons the NMM contribu- 
tion is given by the formula [SI [7] 
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where Ei, is the energy of the incident neutrino, and dis- 
plays a l/T enhancement at low energy transfer. The 



standard electroweak contribution is constant in T at 
> T: 
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In what follows we refer to these two types of contribu- 
tion to the scattering as, respectively, the magnetic and 
the weak. 

The current experiments have reached threshold val- 
ues of T as low as few keV and are likely to further 
improve the sensitivity to low energy deposition in the 
detector. At low energies however one can expect a 
modification of the free-electron formulas ((Ij and ^ 
due to the binding of electrons in the germanium atoms, 
where e.g. the energy of the Ka line, 9.89 keV, indicates 
that at least some of the atomic binding energies are 
comparable to the already relevant to the experiment 
values of T. Thus, a proper treatment of the atomic ef- 
fects in neutrino scattering is necessary and important 
for the analysis of the current and even more of the 
future data with a still lower threshold. For the first 
time this problem was addressed in Ref. [5], where a 
2-3 times enhancement of the electroweak cross section 
in the case of ionization from a Is state of a hydrogen- 
like atom with nuclear charge Z had been numerically 
determined at neutrino energies E^, ^ aZmt? . Subse- 
quent numerical calculations within the Hartree-Fock- 
Dirac method for ionization from inner shells of various 
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atoms showed much lower enhancement 5 — 10%) of 
the electroweak contribution [9l [TOj [11] . The interest to 
the role of atomic effects was renewed in several recent 
papers, which however are ridden by a 'trial and error' 
approach. The early claim [12] of a significant enhance- 
ment of the NMM contribution by the atomic effects has 
been later disproved [131 [13] and it was argued [131 [H] 
that the modification of the formulas ^ and ([2]) by the 
atomic binding effects is insignificant down to very low 
values of T. It has been subsequently pointed out [T5] 
that the analysis of Ref. [T^ is generally invalidated in 
multi-electron systems, including atoms with Z > 1. 
Furthermore, the analysis of Ref. [T3] is also generally 
invalidated by singularities of the relevant correlation 
function in the complex plane of momentum transfer^^ , 
so that the claimed behavior of the cross section at low 
T applies only in the semiclassical limit. 

In this paper we revisit the subject of neutrino scat- 
tering on atoms at low energy transfer. We aim at de- 
scribing this process at T in the range of few keV and 
lower, so that the motion of the electrons is considered 
as strictly nonrelativistic. Also in this range the energy 
of the dominant part of the incident neutrinos from the 
reactor is much larger than T and we thus neglect any 
terms whose relative value is proportional to T / E^. Fur- 
thermore any recoil of the germanium atom as a whole 
results in an energy transfer less than 2E^/MGe, which 
at the typical reactor neutrino energy is well below the 
considered here keV range of the energy transfer. Thus 
we formally set the mass of the atomic nucleus to in- 
finity and neglect any recoil by the atom as a whole. 
In particular, under these conditions the interaction of 
the neutrino with the nucleus can be entirely neglected, 
and only the scattering on the atomic electrons is to be 
considered. 

The kinematics of the scattering of a neutrino on 
atomic electrons is generally characterized by the com- 
ponents of the four- momentum transfer, the energy 
transfer T and the spatial momentum transfer q, from 
the neutrino to the electrons with two rotationally in- 
variant variables being T and q = |q|. At small T the 
electrons can be treated nonrelativistically both in the 
initial and the final state, so that the process is that 
of scattering of an NMM in the electromagnetic field 
A = (^o,A) of the electrons: Ao{q) = v47rap(q)/(j^, 
A(q) = ^/4naj{q) / q"^ , where p{q) and j(q) are the 
Fourier transforms of the electron number density and 
current density operators, respectively. 



^'The flaws in the momentum-transfer dispersion relation and 
sum rules of Ref. |13| are corrected in Ref. |18| . 
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and the sums run over the positions of all the Z elec- 
trons in the atom. 

In this limit the expression for the double differential 
cross section is given by |15j 



dTdq^ 
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where S{T,q^), also known as the dynamical structure 
factor [in], and R{T,q^) are 

^(T,g2) ^^5(T-i?„ + i?o)|(n|p(q)|0)^ (6) 

n 

i?(T, q^) ^Y.^{T^ i?„ + i?o) I (nbl (q)|0) ^ (7) 

n 

with j± being the j component perpendicular to q and 
parallel to the scattering plane, which is formed by 
the incident and final neutrino momenta. The sums 
in Eqs. ([6]) and ([7]) run over all the states \n) with ener- 
gies En of the electron system, with |0) being the initial 
state. 

Clearly, the factors S{T, q^) and i?(T, q^) are related 
to respectively the density-density and current-current 
Green's functions 
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with H being the Hamiltonian for the system of elec- 
trons. For small values of q, in particular, such that 
q ^ T, only the lowest-order non-zero terms of the ex- 
pansion of Eqs. (fTU)) and (fTTj) in powers of q^ are of 
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relevance (the so-called dipole approximation). In this 
case, one has [15] 

R{T,q^)^^SiT,q^). (12) 

q 

Taking into account Eq. p^ . the experimentally 
measured single-differential inclusive cross section is, to 
a good approximation, given by (see e.g. in Refs. [131 
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The standard electroweak contribution to the cross 
section can be similarly expressed in terms of the same 
factor S{T,q^) JJj as 



^^^(l+4sin^. 
dT 47r ^ 
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where the factor S{T, q^) is integrated over q^ with a 
unit weight, rather than as in Eq. ([T^. 

The kinematical limits for q^ in an actual neutrino 
scattering are explicitly indicated in Eqs. (|T3)) and (IT4l) . 
At large E^, typical for the reactor neutrinos, the upper 
limit can in fact be extended to infinity, since in the dis- 
cussed here nonrelativistic limit the range of momenta 
~ E,^ is indistinguishable from infinity. The lower limit 
can be shifted to = 0, since the contribution of the re- 
gion of q^ < T^ can be expressed in terms of the photo- 
electric cross section [TB] and is negligibly small (at the 
level of below one percent in the considered range of T) . 
For this reason we henceforth discuss the momentum- 
transfer integrals in Eqs. (jl3p and (jl4p running from 



to q^ 



h{T)^ rS{T,q')^, 
Jo ^ 
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For a free electron, which is initially at rest, the 
density-density correlator is the free particle Green's 
function 

so that the dynamical structure factor is given by 
S{FE){T,q'^) — S{T - q^/2m), and the discussed here 
integrals are in the free-electron limit as follows: 

^d(f_ _ l_ 



(FE) 



S(FE){T,q')dq^ ^2m.. (17b) 



It is readily seen that these expressions, when used in 
the formulas (jl3p and (fT4)) . result in the free-electron 
cross section in Eqs. (H]) and 

Let us consider the scattering on just one bound elec- 
tron. The Hamiltonian for the electron has the form 
H ~ /2ra + V{r), and the density-density Green's 
function from Eq. ([S]) can be written as 



F(T,q2) = (0 e-*i'"[T-i7 + £;o]"^e*'i"" 



T-f-^-H+E:;" 

2m m 



, 



(18) 



where the infinitesimal shift T — > T — ie is implied. 

Clearly, a nontrivial behavior of the latter expression 
in Eq. (IT51) is generated by the presence of the opera- 
tor pq/m in the denominator, and the fact that it does 
not commute with the Hamiltonian H . Thus an ana- 
lytical calculation of the Green's function as well as the 
dynamical structure factor and the momentum-transfer 
integrals is feasible in only few specific problems. In 
particular, such a calculation for ionization from the Is, 
2s, and 2p hydrogen- like states shows that the deviation 
of the discussed integrals (|15p from their free values are 
very small [18]: the largest deviation is exactly at the 
ionization threshold, where, for instance, each of the Is 
integrals is equal to the free-electron value multiplied 
by the factor (1 — 7e^'*/3) ~ 0.957, and in the 2s and 
2p cases the departure from the free-electron behavior 
is even smaller. 

The problem of calculating the integrals (ITSt how- 
ever can be solved in the semiclassical limit, where one 
can neglect the noncommutativity of the momentum 
p with the Hamiltonian, and rather treat this opera- 
tor as a number vector. Taking also into account that 
(i/ — £^o)|0) = 0, one can then readily average the latter 
expression in Eq. (IT51) over the directions of q and find 
the formula for the dynamical structure factor: 



S{T,q^ 
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pq 
m 



(19) 



where p = |p| and 9 is the standard Heaviside step 
function. The expression in Eq. (jl9|) is nonzero only 
in the range of q satisfying the condition ~pq/m < 
(17a) T — q'^ /2m < pq/m, i.e. between the (positive) roots 
of the binomials in the arguments of the step functions: 
q.m^n = ^/2mT + p'^ - p and qmax = \/ '^mT + p^ + p. 
One can notice that the previously mentioned 'spread 
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and shift' of the peak in the dynamical structure func- 
tion in this hmit corresponds to a flat pedestal between 
qmin and qmax- The calculation of the integrals (jlSp 
with the expression (1191) is straightforward, and yields 
the free-electron expressions (jl7l) for the discussed here 
integrals in the semiclassical (WKB) limit: 

/j'^^^) = i = 2m. (20) 

The difference from the pure free-electron case however 
is in the range of the energy transfer T. Namely, the 
expressions (^0)) are applicable in this case only above 
the ionization threshold, i.e. at T > |i5o|. Below the 
threshold the electron becomes 'inactive'. 

It is instructive to point out that the validity of the 
result in Eq. (PH)) is based on the semiclassical approx- 
imation and is not directly related to the value of the 
energy T. In particular, for a Coulomb interaction the 
WKB approximation is applicable at energy near the 
threshold |17]. For T exactly at the threshold, T = -Ea, 
the criterion for applicability of the semiclassical ap- 
proach in terms of the force F = |F| ~ |[p._ff]| acting 
on the electron and the momentum p of the electron is 
that [17] the ratio of the characteristic values mF/p^ 
is small. For the excitation of a state with the prin- 
cipal number n one has \F\ = a/r^ ~ m^a^n"* and 
p ^ ma/n, so that m\F\/p'^ ^ 1/n. Thus the applica- 
bility of a semiclassical treatment of the ionization near 
the threshold improves for initial states with large n. As 
previously mentioned, the modification of the integrals 
(fTS]) by the binding is already less than 5% for n = 1, 
so that we fully expect this deviation to be smaller for 
the higher states, and even smaller at larger values of 
T above the threshold due to the approach to the free- 
electron behavior a,t T ^ Eq. 

We believe that the latter conclusion explains the 
so-called stepping behavior observed empirically [10^ in 
the results of numerical calculations. Namely the calcu- 
lated cross section da/dT for ionization of an electron 
from an atomic orbital follows the free-electron depen- 
dence on T all the way down to the threshold for the 
corresponding orbital with a very small, at most a few 
percent, deviation. This observation led the authors of 
Ref. [TOj to suggest the stepping approximation for the 
ratio of the atomic cross section (per target electron) to 
the free-electron one: 

«^>-(5^4i:".«(^-iE.i). (21) 

where the sum runs over the atomic orbitals with the 
binding energies Ei and the filling numbers n;. Clearly, 
the factor /(T) simply counts the fraction of 'active' 



electrons at the energy T, i.e. those for which the ion- 
ization is kinematically possible. For this reason we refer 
to f{t) as an activation factor. 

In considering the neutrino scattering on actual 
atoms one needs to evaluate the dependence of the num- 
ber of active electrons on T. The energies of the inner 
K, L, and M orbitals in the germanium atom are well 
known (see e.g. Ref. [TT] and references therein) and 
provide the necessary data for a description of the neu- 
trino scattering by the stepping formula (j2ip down to 
the values of the energy transfer T in the range of the 
binding of the M electrons, i.e. at T > \Em\ ^ 0.1 keV. 
The corresponding steps in the activation factor are 
shown in Fig. [TJ 
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Fig. 1. The activation factor / for germanium in the 
stepping approximation with the actual energies of the 
orbitals (solid line) and its value in the Thomas-Fermi 
model (dashed). 



Let us show that the stepping approximation ([?!]) 
can also be explained in terms of quantum corrections to 
the activation factor given by the Thomas- Fermi atomic 
model (see e.g. Ref. [ITJ). The latter describes electrons 
in an atom as a degenerate free electron gas in a master 
potential 4)(r) filling the momentum space up to the zero 
Fermi energy, i.e. up to the momentum po{r) such that 
p^/2m — ecj) = 0. The electron density n{r) — pg/(37r^) 
then determines the potential (j){r) from the usual Pois- 
son's equation, thus resulting in a self-consistent proce- 
dure. In the discussed picture at an energy transfer T 
the ionization is possible only for the electrons whose en- 
ergies in the potential are above — T, i.e. with momenta 
above prir) with p"^/ 2m — ecj) — ~T. The electrons with 
lower energy are inactive. Calculating the density of the 
inactive electrons as pij./(37r^) and subtracting their to- 
tal number from Z , one readily arrives at the formula 
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for the activation factor, i.e. the effective fraction of the 
active electrons Z^s/Z as a function of T: 



f{T) = 



Zes{T) 



= 1 - 





'x{x) 


T' 


Jo 


X 





3/2 



(22) 

where x(a;) is the Thomas-Fermi function, well 
known and tabulated, of the scaling variable 
X = 2(4/37r)^/'^TOQ;Z^/'^, the energy scale Tq is 
given by 



4 \ 



2/3 



Tn = 2 ( — j ma^Z^/^ w SO.SZ''/^ eV, 



(23) 



and, finally, X(){T) is the point where the integrand be- 
comes zero, i.e. corresponding to the radius beyond 
which all the electrons are active at the given energy T. 
The energy scale Tq in germanium (Z=32) evaluates to 
To ss 3.1 keV. The activation factor for germanium cal- 
culated from the formula ([^^ is shown by the dashed 
line in the plot of Fig. 2. One can see that the stepping 
activation factor (j2ip mimics upon average over the en- 
ergy intervals between the electron shells in germanium 
the Thomas-Fermi result. Thus, it can be considered as 
refinement of the latter due to account for the quanti- 
zation of the electron binding energies. 

It should be remarked that the discussed statistical 
model is known to approximate the average bulk prop- 
erties of the atomic electrons with a relative accuracy 
0{Z~^/'^) and as long as the essential distances r sat- 
isfy the condition Z^^ ^ mar ^ 1, which condition in 
terms of the scaling variable x reads as Z~^^^ ^ x <C 
Z^^^. In terms of the formula for the number of 
active electrons, the lower bound on the applicability 
of the model is formally broken at T ~ Z^^'^Tq, i.e. at 
the energy scale of the inner atomic shells. However the 
effect of the deactivation of the inner electrons is small, 
of order Z~^ in comparison with the total number Z of 
the electrons. On the other hand, at low T, including 
the most interesting region of T ^ Tq, the integral in 
Eq. (I22p is determined by the range of x of order one, 
where the model treatment is reasonably justified. 

We have considered the scattering of neutrinos on 
electrons bound in atoms. Our main finding is that the 
differential over the energy transfer cross section given 
by the free-electron formulas ([T]) and (0) and the step- 
ping behavior of the activation factor given by Eq. (pij) 
provides a very accurate description of the neutrino- 
impact ionization of a complex atom, such as germa- 
nium, down to quite low energy transfer. The devia- 
tion from this approximation due to the onset of the 
ionization near the threshold is less than 5% (of the 
height of the step) for the K electrons, if one applies 



the analytical behavior of this onset that we find for the 
ground state of a hydrogen-like ion. We also find that 
the free-electron expressions for the cross section are not 
affected by the atomic binding effects in the semiclassi- 
cal limit. For this reason we expect that the deviation 
of the actual onset from a step function at the threshold 
for ionization of higher atomic orbitals is even smaller 
than for the ground state, since the motion in the higher 
states is closer to the semiclassical limit. Thus, our an- 
alytical results explain the numerically determined be- 
haviors of the electroweak and magnetic contributions 
to the neutrino-impact ionization of various atomic tar- 
gets [iiiniin]. 
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